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Abstract 

Background: Field observations and a few physiological studies have demonstrated that peach embryogenesis 
and fruit development are tightly coupled. In fact, attempts to stimulate parthenocarpic fruit development by 
means of external tools have failed. Moreover, physiological disturbances during early embryo development lead to 
seed abortion and fruitlet abscission. Later in embryo development, the interactions between seed and fruit 
development become less strict. As there is limited genetic and molecular information about seed-pericarp cross- 
talk and development in peach, a massive gene approach based on the use of the uPEACH 1.0 array platform and 
quantitative real time RT-PCR (qRT-PCR) was used to study this process. 

Results: A comparative analysis of the transcription profiles conducted in seed and mesocarp (cv Fantasia) 
throughout different developmental stages (S1, S2, S3 and S4) evidenced that 455 genes are differentially 
expressed in seed and fruit. Among differentially expressed genes some were validated as markers in two 
subsequent years and in three different genotypes. Seed markers were a LTP1 (lipid transfer protein), a PR 
(pathogenesis-related) protein, a prunin and LEA (Late Embryogenesis Abundant) protein, for SI, S2, S3 and S4, 
respectively. Mesocarp markers were a RD22-like protein, a serin-carboxypeptidase, a senescence related protein 
and an Aux/IAA, for S1, S2, S3 and S4, respectively. 

The microarray data, analyzed by using the HORMONOMETER platform, allowed the identification of hormone- 
responsive genes, some of them putatively involved in seed-pericarp crosstalk. Results indicated that auxin, 
cytokinins, and gibberellins are good candidates, acting either directly (auxin) or indirectly as signals during early 
development, when the cross-talk is more active and vital for fruit set, whereas abscisic acid and ethylene may be 
involved later on. 

Conclusions: In this research, genes were identified marking different phases of seed and mesocarp development. 
The selected genes behaved as good seed markers, while for mesocarp their reliability appeared to be dependent 
upon developmental and ripening traits. Regarding the cross-talk between seed and pericarp, possible candidate 
signals were identified among hormones. 

Further investigations relying upon the availability of whole genome platforms will allow the enrichment of a 
marker genes repertoire and the elucidation of players other than hormones that are involved in seed-pericarp 
cross-talk (i.e. hormone peptides and microRNAs). 
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Background 

Peach fruit development is tightly connected to embryo- 
genesis. Fruit growth displays a double sigmoid pattern 
in which four stages named SI, S2, S3 and S4 can be 
distinguished [1]. The early part of SI is characterized 
by cell division and enlargement lasting about two 
weeks, followed by cell enlargement. The slowdown in 
growth that occurs at S1/S2 transition is followed by 
endocarp lignification (pit hardening), which lasts for 
12-15 days from the middle of S2 to its end. S3 starts 
with a resumption of growth mainly due to cell enlarge- 
ment, thus generating the second exponential phase. 
Maturation is completed by the end of S3 and followed 
by ripening (S4). The four fruit developmental phases 
are determined using a mathematical model based on 
first derivative of the growth curve [1]. Identification of 
the growth phases is important both for developmental 
studies and for precision farming. However, the only 
easily detectable event is the end of pit hardening mark- 
ing the S2/S3 transition, because the phase length is 
affected by both genotype (early, middle and late ripen- 
ing varieties) and environmental cues. A continuous 
growth model reassessment is therefore required. 
Accordingly, the identification of developmental phase 
organ-specific molecular markers would be of great 
importance for scientific and practical purposes. 

Seed development, necessary for fruit set [2], is char- 
acterized by a fast endosperm growth that starts imme- 
diately after fertilization concurrently with the nucellus 
re-absorption, and lasts until the beginning of endocarp 
lignification, when the seed reaches its final size. At the 
end of pit hardening, seed volume is mainly made up of 
endosperm and the embryo is at the heart stage. There- 
after, embryo growth resumes and cotyledon develop- 
ment is paralleled by endosperm re-absorption. Seed 
maturation is characterized by lipids accumulations [3], 
synthesis of specific late embryogenesis abundant (LEA) 
proteins and dehydration. Attempts to stimulate parthe- 
nocarpic fruit development by hormone applications 
resulted as being ineffective. Moreover, seed abnormal- 
ities at the early stages of development (SI and S1/S2 
transition stages) lead to abortion and fruitlet abscission 
[4]. Later, (late S2, S3 and S4), the relationships between 
fruit development and embryogenesis become less strict. 
This is the case for early ripening varieties characterized 
by the uncoupling of fruit development and embryogen- 
esis. In fact, at harvest, seed development is still in pro- 
gress and a long way from maturity. Seed presence is 
always necessary to achieve normal fruit development 
even if embryo development is incomplete [5], Apart 
from the above observations, molecular-genetic informa- 
tion on the relationship between fruit and seed develop- 
ment is scarce. Cross-talk between the two organs may 
involve different components of the signaling network, 



such as hormones, transcription factors (TFs) and other 
signaling molecules, playing either direct or indirect 
roles. 

Concerning hormones, parthenocarpic fruit develop- 
ment in some species is induced by applications of 
auxin or cytokinins (CKs), or gibberellins (GAs), or hor- 
mone blends [6]. Molecular approaches have confirmed 
the role played by hormones, especially auxins [7], 
Investigations in Arabidopsis identified a mutant, named 
fwf (fruit without fertilization), with a normal silique 
development even in the absence of seeds [8]. Double 
mutant analysis (fwf gal-4, fwf gai, fwf spy, fwf ats) 
pointed out that FWF negatively affected GA biosynth- 
esis and GA and auxin signal transduction. The FWF 
protein may interact with TFs such as Fruitful (FUL) 
and Aberrant Testa Shape (ATS), members of the 
MADS-box family, and Scarecrow -SCR- type, which are 
all involved in cell division [8]. Additional TFs have 
been identified, some of which are related to hormone 
action, actively transcribed along peach fruit develop- 
ment and ripening ([9,10]). Orthologues of these TFs 
are also expressed in true (silique and berry) and false 
(pome and strawberry) fruits, supporting the hypothesis 
that different fruit types share common regulatory ele- 
ments [11]. High throughput analysis conducted in Ara- 
bidopsis showed that some TFs are shared by seed and 
fruit [12]. 

Taking this information into account, peach seed and 
fruit transcriptomes were explored throughout develop- 
ment by means of a massive gene approach based on 
the use of the uPEACH 1.0 array platform and quantita- 
tive real time RT-PCR (qRT-PCR). The research identi- 
fied genes marking organ/tissue developmental phases, 
as well as candidate signals (hormones and TFs) that 
may trigger the cross-talk between fruit and seed. 

Results 

Seed and fruit growth pattern 

Fruit growth analysis was performed on cv Fantasia and 
assumed as a reference (Figure 1). In this genotype fruit 
development and ripening are completed in 135-140 
days after full bloom (DAFB). Growth dynamics display 
the typical double sigmoidal pattern in which four 
developmental stages have been identified according to 
the first derivative. SI, S2, S3 and S4 lasted for 45, 32, 
33 and 17 days, respectively. Pit hardening (PH) started 
60 DAFB and was completed by the S2/S3 transition. 
The seed derives from the fertilized ovule and the initial 
increase in length (Figure 1) is due to the rapid nuclear 
division of the endosperm responsible for embryo sac 
expansion. Endosperm cellularization starts 40 DAFB 
and is completed by the beginning of PH. The embryo 
develops very slowly in the early stages (SI and S2), 
reaching a length of about 40-60 um. Later, at the S2/S3 
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Figure 1 Fruit and seed growth pattern (cv Fantasia) Fruit 
growth (red) is expressed as cross diameter while length is used for 
seed (blue) and embryo (green) development. Difference in length 
between seed and embryo represents endosperm, integuments and 
nucellus being a minimal part of the seed. Fruit developmental 
cycle has been divided into 4 main stages (S1 to S4) according to 
the first derivative of the fruit growth curve. The yellow horizontal 
ine indicates pit hardening. Sampling dates are marked by black 
arrows. The simple loop microarray experimental design is outlined 
on the right. For the microarray expression analyses, seed (S) and 
mesocarp (M) tissues at S1 and S2I, and S3 and S4 were pooled, 
and defined as early (ES and EM) and late (LS and LM) 
development, respectively. The comparison has been made 
between different developmental stages (LS/ES and LM/EM) within 
the organs and between the two organs (ES/EM and LS/LM) within 
the developmental stage. 



transition, it resumes development reaching its final size 
by the middle of S3. The morphological completion of 
development is followed by maturation and desiccation. 

Identification of marker genes 

RNAs extracted before (E, early development) and after 
(L, late development) pit hardening have been used for 
microarray transcriptome analyses in order to identify 
genes possibly involved in seed-pericarp cross-talk or 
useful as organ and developmental phase molecular 
markers. Data obtained from the microarray analyses 
were handled either as single comparisons, i.e. late seed 
vs. early seed (LS/ES), late mesocarp vs. early mesocarp 
(LM/EM), within each hybridization or by combining 
the whole set of data, thus also including ES/EM and 
LS/LM (see Figure 1 insert). The microarray expression 
data (see Additional file 1), validated by means of qRT- 
PCR on 29 randomly selected genes, showed a Pearson 
correlation coefficient ranging, in the four comparisons, 
from 0.79 to 0.84 (see Additional file 2). 

With the single comparison analyses, among the 360 
differentially expressed genes within the two organs at 
early and late development (Figure 2A), 174 and 151 
were differentially expressed only in seed (groups A and 
B) and mesocarp (groups C and D), respectively. Of the 




B 




Figure 2 Genes differentially expressed according to the 
developmental stage of the organ. Venn diagrams were used to 
visualize genes differentially expressed in the microarray 
experiments. Comparisons between early (E) and late (L) 
development (panel A), and seed (S) and mesocarp (M) (panel B), 
were made by means of a direct comparison approach (LS/ES and 
LM/EM in A; ES/EM and LS/LM in B). Arrowhead orientation 
indicates up ( A ) and down (▼) regulation. The letters inside the 
sectors are tags for the identification of the genes listed in 
Additional file 1 . 



seed differentially expressed genes, 108 and 66 were 
more transcribed at early (group B) and late develop- 
ment (group A), respectively. Four genes, shared by seed 
and mesocarp, were more actively transcribed at late 
development (group E), while an additional four showed 
the opposite trend of expression, being induced in LS 
and repressed in LM (group H). In addition to the 108 
genes more abundantly transcribed in ES (group B), 22 
were also expressed in EM (group G), while 5 were 
abundant in ES and EM (group F). Among the meso- 
carp differentially expressed genes, 101 and 50 were 
more transcribed in EM (group D) and LM (group C), 
respectively. Taking the comparison between seed and 
mesocarp (ES/EM and LS/LM) into account, 341 genes 
were differentially expressed in the two organs (Figure 
2B). Among these, 133 and 151 were differentially 
expressed only at early (groups I and L) and late (groups 
M and N) development, respectively. Considering the 
differentially-expressed genes at early development, 40 
mRNAs were more abundant in seed (group I) and 93 
in mesocarp (group L), while among the late develop- 
ment ones, 97 were more abundant in seed (group M) 
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and 54 in mesocarp (group N). Of the 57 remaining 
genes, 17 and 35 transcripts were always more abundant 
in seed (group O) and mesocarp (group Q), respectively, 
and 5 displayed an opposite pattern, being more (3) or 
less (2) abundant in ES (group R) and ML (group P). 
Annotations of genes included in Figure 2 are reported 
with microarray expression data in Additional file 1. 

Based on the above microarray analysis, putative mar- 
kers were searched to find those that meet the following 
criteria: a) moderately to highly expressed in only one 
organ (seed or mesocarp), b) highly expressed/not 
expressed at specific developmental stage/s (SI to S4). 
According to these criteria, 50 potential marker genes, 
chosen among those differentially expressed in the 
microarray, were selected and tested by means of qRT- 
PCR in leaf, flower (data not shown), seed and mesocarp 
at five developmental stages in cv Fantasia (Figure 3). 
These detailed expression profiles allowed the identifica- 
tion of eight genes best fulfilling the ideal marker 



criteria. For seed development, ctg3431, coding for a 
lipid transfer protein (LTP), ctgl026, coding for a patho- 
genesis related (PR) protein, ctgl540, coding for a pru- 
nin, and ctg3563, coding for a late embryogenesis 
abundant (LEA) protein, have been chosen as SI, S2, S3 
and S4 markers, respectively. Concerning mesocarp 
development, ctg2909, coding for a RD22-like protein, 
ctgl751, coding for serine carboxypeptidase, ctgl823, 
coding for a senescent associated protein, and ctg57, 
coding for an AUX/IAA protein, have been selected as 
SI, S2, S3 and S4 markers, respectively (Figure 3). The 
function as stage markers has been confirmed on the 
same genotype for an additional growing season (Addi- 
tional file 3). 

A further validation of the selected genes was per- 
formed in two additional genotypes (cv Springcrest and 
the slow ripening - sir - selection) differing for the 
dynamics of seed and fruit development. In Springcrest, 
fruit ripening occurred after 86 DAFB (Figure 4A), 




ctg name annotation 



1357 MADS-box TF 
2195 MADS-box TF 
3430 LTP 



3431 LTP 




1026 PR-4B 




1151 ACS3 




2116 ERF2 




2916 SSDH 




1434 RuBisCO 




1543 prunin 




1540 prunin 




544 GA2 oxidase 



4056 IAA amidohydrolase 
5138 Ala glyoxylate aminotransfera 
3563 LEA 
2025 ETR2 

2370 hist phosp protein 
2499 EREBP 
2655 SAUR 
4163 unknown 
2909 RD-22 
1586 unknown 
974 RD-22 
349 PIP 
2045 senescence protein 

945 CAB | 75 

2568 ACS2 | 100 

1024 catalase 
2605 calmodulin 

1257 invertase/pectin methylesterase inhibitor 
2902 sorbitol transporter 
1751 serine carboxypeptidase 
1421 WRKYTF 
1436 ETR1 

56 ABA WDS superfamily 
1528 dehydrin 
1819 senescence protein 
1823 senescence protein 
2599 membrane protein 
1391 GAST-like protein 

356 ERS 
298 bZIP 
358 Aux/IAA 
64 AC01 

57 Aux/IAA 
1993 GH3 

1112 unknown DUF538 
939 EXP1 
941 EXP2 
489 ACS1 

Figure 3 Selection of developmental stage and organ specific marker genes Identification of putative marker genes was performed by 
selecting some of those differentially expressed in the microarray analyses and further validated by means of qRT-PCR. This detailed expression 
profiling allowed the selection of those genes that best fitted the ideal marker characteristics as indicated in the Methods section. Expression 
profiles of 50 genes were measured in seed and mesocarp at five different developmental stages (S1 to S4). Expression values, as indicated in 
the insert, are related to the highest expression of each gene (100% blue). Genes have been manually ordered according to their expression 
profiles. Grey shading highlights genes selected as markers. 
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Figure 4 Dynamics of fruit and seed growth in Fantasia, Springcrest and sir A) Fruit growth curves are expressed as cross diameter (mm) 
for Fantasia (the reference genotype; red triangles), Springcrest (the early ripening genotype; blue squares) and sir (the slow ripening genotype; 
green circles). In the lower part of the panel, the arrowheads indicate the timing of sampling for the 3 cvs and the developmental stage is 
indicated within each arrow. B) Dynamics of seed development in Springcrest (left) and Fantasia (right) related to the fruit developmental stages. 
Seed development in sir is similar to that reported for Fantasia. Relative abundance of nucellus, integuments and endosperm (blue) and embryo 
(red) points out that in Springcrest, at fruit harvest, embryo development is a long way from maturity, while in sir, in spite of the block of fruit 
ripening, the completion of embryo development parallels that of Fantasia and the seed is viable. 



when seed development was still in progress (Figure 4B). 
At the end of the growing season (taking cv Fantasia as 
a reference), sir showed a fully developed seed (Figures 
4A and 4B), while the mesocarp development was 
blocked at stage S3. 

As regards seed markers, ctg3431, coding for a LTP, 
clearly marked the SI stage for both Fantasia and sir, 
while in Springcrest its expression decreased only at S3 
stage (Figure 5A). A PR protein encoding gene, ctgl026, 
has been selected for the S2 stage. The highest expres- 
sion level was found in the seed of cv Fantasia, peaking 
at early S2 and decreasing thereafter, as in Springcrest. 
In sir, its expression was broader, being relevant also at 
SI and S2II (87 and 86% of S2I, set as 100%, respec- 
tively; Figure 5B). A prunin, the main seed storage 



protein in Prunus spp., encoded by ctgl540, is a good 
marker for S3 seed development only in Fantasia. In 
fact, different amounts and kinetics of its transcript 
accumulation were observed in the other two genotypes. 
In Fantasia, accumulation started between S2I and S2II 
and increased up to a maximum at S3, decreasing there- 
after, whereas in sir and Springcrest transcript accumu- 
lation was delayed, becoming detectable at S3 in the 
former and S4 in the latter (Figure 5C). The expression 
of the gene encoding a LEA protein (ctg3563) became 
detectable at S2II in Fantasia and peaked at S4. A simi- 
lar pattern was observed in sir, although the transcript 
only started to be detectable at S3. In Springcrest, it was 
detectable only at S4, at levels lower than in the other 
two genotypes (Figure 5D). The level of expression of 
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Figure 5 Validation of developmental stage and organ specific 
markers in mesocarp and seed of three genotypes. Expression 
pattern, assessed by qRT-PCR, of seed (dashed lines) and mesocarp 
(solid lines) molecular markers of Fantasia (red triangles), Springcrest 
(blue squares) and sir (green circles), at five developmental stages 
(S1 to S4). Transcript levels are measured as means of normalized 
expression ± SEM of three technical replicates. 



the four genes in mesocarp was very low throughout 
development and comparable in the three cvs (Figures 
5A-D). 

As regards mesocarp, ctg2909, coding for an RD22- 
like protein, had maximum expression at SI and early 
S2 (i.e. S2I, Figure 5E). In Fantasia and sir the expres- 
sion decreased already at S2II (28% and 32% of the max- 
imum in Fantasia and sir, respectively), while in 
Springcrest its expression was still high (96%) at S2II. 

A serine carboxypeptidase (ctgl751) was chosen as a 
marker for the S2 developmental stage. In Fantasia, the 
transcript was undetectable at SI, at basal level at S2I, 
peaked sharply at S2II, and then declined at S3 and S4. 
Also in the other two varieties the transcript was unde- 
tectable at SI, but its expression, already high at S2I, 
slightly increased at S2II and remained at high levels at 
S3, decreasing at S4 (Figure 5F). The expression of 
ctgl823, encoding a senescence related protein, had a 
maximum in Fantasia at S3 (100%), while expression 
levels were much lower in the previous and following 
stages (29 and 9% at S2II and S4, respectively). Although 
its expression was relatively high (50%) at S2I, it may be 



considered a good S3 marker. In Springcrest, the expres- 
sion was generally low at all stages, with a maximum at 
S2II. In the sir genotype, the accumulation of ctgl823 
transcripts steadily increased during the early phases up 
to a maximum at S2II. Although slightly decreasing 
thereafter, the ctgl823 mRNA was also abundant at S3 
and S4 (60 and 74% of S2II, respectively) (Figure 5G). 
S4 stage is clearly identified by the expression of ctg57, 
coding for an Aux/IAA protein. In Fantasia, the expres- 
sion at S3 is about 6% of that measured at S4 and 
almost undetectable in early phases. In Springcrest its 
expression is also almost undetectable at SI, S2I and 
S2II, but at S3 it is already half of that measured at S4. 
In sir, although maximum expression is at S4, the tran- 
scripts accumulated at very low levels (5% of Fantasia) 
(Figure 5H). 

In agreement with their being mesocarp markers, all 
the selected genes are almost undetectable in seed (Fig- 
ure 5E-H) with the exception of ctgl823 in sir. 

Hormones and TFs in seed fruit cross-talk 

Hormone-related genes possibly involved in cross-talk 
between the two organs were identified among those 
spotted on the microarray based upon the list of hormo- 
nal indexes available for Arabidopsis ([13]; TAIR web- 
site). The portion of hormone responsive genes in 
Arabidopsis ranges between 3.8 and 9.4% of the whole 
transcriptome (TAIR 10 vers., 27,416 genes), depending 
on the hormone considered (Additional file 4). For 
uPEACH 1.0 (4,806 targets), the portion of hormone 
responsive genes parallels that of Arabidopsis, ranging 
from 3.8 to 9.8% with values for each hormone class 
comparable to those calculated for Arabidopsis. An irre- 
levant bias may therefore be assumed to exist when 
peach expression data are used for HORMONOMETER 
analysis [13]. In addition, it could be assumed that the 
same proportion might be expected if a whole genome 
array were used. 

A heat map was produced by considering the follow- 
ing subsets of genes for each hormone (Figure 6): i) 
genes involved in signal transduction (ST), ii) hormone- 
responsive genes (H), iii) genes with hormone-specific 
responsiveness (SRG), iv) hormone-responsive genes 
encoding TFs (TFs), and v) genes encoding TFs with 
hormone-specific responsiveness (sTFs). The subset i) 
was identified using the classification of Arabidopsis 
orthologs obtained from TAIR GO terms and AHD 
classification lists (available at http://ahd.cbi.pku.edu.cn/; 
[14]), and was then analyzed by averaging the log ratios, 
while the other subsets were used for the HORMON- 
OMETER analyses [13]. 

Concerning auxin and intra-organ comparisons (LS/ES 
and LM/EM), a weak activation of ST was observed in 
LS with respect to ES, paralleled by a partial correlation 
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Figure 6 Heat map showing the relationship between the expression of signal transduction and hormone target genes Panel A The 

heat map was produced by considering the genes involved in the signal transduction (ST) for auxin (AUX), cytokinin (CK), gibberellic acid (GA), 
abscissic acid (ABA) and ethylene (C 2 H 4 ). HORMONOMETER data were grouped into hormone-responsive genes (H), genes with hormone-specific 
responsiveness (SRG), hormone-responsive genes encoding TFs (TFs), and genes encoding TFs with hormone-specific responsiveness (sTFs). For 
each hormone, the following comparisons have been analyzed: LS/SE, LM/EM, ES/EM and LS/LM. Panel B Color codes for ST genes and 
hormone-responsive genes (HORMONOMETER). For ST, red and green represent up- and down-regulation, respectively. In the HORMONOMETER, 
orange (value = 1), white (value = 0), and blue (value = -1) indicate a complete correlation, no correlation, or anti-correlation, respectively, in 
terms of direction and intensity of the hormone index with the queried experiment [13]. 



with the overall reference hormone indexes, whereas a 
partial anti-correlation was observed when auxin-specific 
hormone indexes, TF- and specific TF-encoding targets 
were used. In the mesocarp, a marked up-regulation of 
ST subset was evidenced in LM, and a good correlation 
was shown in the same sample both considering the 
overall hormone indexes and all the other gene subsets. 
As regards inter-organ comparisons, a decreased tran- 
scription of ST elements was always observed in the seed, 
paralleled by an anti-correlation with the overall hor- 
mone indexes at both early (ES/EM) and late (LS/LM) 
development. However, considering the specific subset, a 
slight correlation was found in the former comparison, 
whereas all the results in the latter one were consistent 
with the overall HORMONOMETER data. 

The intra-organ comparison LS/ES indicated a down- 
regulation of cytokinin (CK) ST elements at late seed 
development, paralleled by an anti-correlation with both 
the overall and specific hormone indexes. However, a 
slight correlation was observed in terms of specific TFs, 
while all TFs appeared not correlated. Concerning the 
mesocarp, a lower activation of ST elements in LM than 
EM was counteracted by a strong correlation with CK 
indexes. CK-specific genes appeared not correlated, 
whereas TFs showed a slight correlation, becoming 
stronger when only the CK-specific TFs were consid- 
ered. As regards inter-organ comparisons, a low activa- 
tion of the signal transduction in ES was counteracted 
by a strong correlation with overall hormone indexes. 
When the analysis was performed with the other sub- 
sets, a significant anti-correlation was observed. Finally, 
during late seed development, despite the higher activa- 
tion of ST elements compared to the mesocarp, a 



general anti-correlation was shown, with the exception 
of specific TFs, which appeared not correlated. 

Considering the gibberellins (GAs)-related expression 
data, the LS/ES comparison demonstrated a good con- 
sistency in signal transduction, and anti-correlation with 
overall and specific transcriptional indexes, and TFs, 
except for the GA-specific TFs, that were not correlated. 
The mesocarp profile was similar except when all TFs 
were considered, the latter analysis showing a robust 
correlation. In the ES/EM inter-organ comparison, a 
depression of the ST pathway in the seed was evidenced. 
The overall HORMONOMETER analysis showed no 
correlation with GA hormone indexes, whereas an anti- 
correlation resulted from the analysis of hormone-speci- 
fic targets. When all the TFs underwent the HORMON- 
OMETER analysis, a strong correlation was shown, 
while specific TFs were not correlated. The most signifi- 
cant data pointed out by the LS/LM comparison con- 
cerned the analysis of GA-specific indexes, showing a 
slight correlation. 

As regards abscisic acid (ABA) and intra-organ com- 
parisons, in spite of a down-regulation of its ST pathway 
during late seed development, a correlation was 
observed in terms of both overall and ABA-specific 
indexes. TFs were basically anti-correlated and not cor- 
related, when considered either as a whole or just the 
specific ones, respectively. In the mesocarp, despite a 
weak up-regulation of the ST elements found in LM, 
there was no significant correlation in any of the HOR- 
MONOMETER analyses. Moving to inter-organ com- 
parison ES/EM, the down-regulation of signal 
transduction occurring in ES paralleled an anti-correla- 
tion found in all the gene sets. In the LS/LM 
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comparison, similar results were obtained in terms of 
both signal transduction and HORMONOMETER. 

Concerning ethylene, no variation was observed 
between LS and ES in terms of expression of genes 
encoding ST elements. In spite of this, a slight correla- 
tion was pointed out by both overall and ethylene-speci- 
fic gene targets. Moreover, TFs were not correlated, 
while specific TFs were slightly anti-correlated. With the 
LM/EM comparison, the hormone signaling pathway 
was up-regulated in LM, paralleled by a partial correla- 
tion of TFs. On the other hand, both the hormone spe- 
cific subsets showed an anti-correlation, stronger in the 
case of TFs. Both inter-organ comparisons (ES/EM and 
LS/LM) displayed a down-regulation of the ST pathway 
in the seed. The HORMONOMETER analyses showed 
no correlation when all targets and all TFs were consid- 
ered, and anti-correlation concerning the specific targets 
and TFs, stronger for the former. Both signal transduc- 
tion and HORMONOMETER results related to jasmo- 
nates, salicylic acid, and brassinosteroids are presented 
and discussed in Additional file 5. 

Discussion 

This research was mainly focused on the relationship 
between seed and pericarp throughout development, 
using a mass gene approach by means of the 
uPEACHl.O [9]. Although this platform was developed 
mainly from late development mesocarp cDNAs, hybri- 
dization analyses and differential expression profiles 
assessed for both early developing mesocarp and seed 
indicate that uPEACHl.O is also a reliable tool for these 
transcriptomic investigations. 

Concerning marker genes, morphological observations 
pointed out that the dynamics of seed development in 
different genotypes is quite synchronous, whereas a 
wide variation exists in the pericarp, affecting not only 
the length of the developmental phases but also impor- 
tant traits related to fruit quality, such as the degree of 
endocarp lignification (cartilaginous endocarp), flesh 
texture (melting/non-melting), sugar/acid ratio, etc. 
Accordingly, the singling out of marker genes specific 
for the same developmental stage is not always unequi- 
vocal for all three studied genotypes. Moreover, since 
seed sampling was referred to the fruit developmental 
stages (SI, S2, S3 and S4), expression data should be 
read taking into account the uncoupling that exists 
between seed and fruit development in Springcrest, an 
early ripening cultivar. 

The ctg3431, marking SI in the seed, encodes a lipid 
transfer protein similar to Arabidopsis LTP1 [15]. Its 
gene expression profile in peach is consistent with Ara- 
bidopsis data, the latter showing that LTP1, along with 
LTP3, LTP4 and LTP6, is expressed at high levels during 
early seed development [16]. The function of this gene 



as an SI marker was confirmed in all the genotypes. 
The delayed decay of transcript accumulation assessed 
in the seed of cv Springcrest has, in fact, to be related 
to the acceleration of mesocarp development in this 
genotype (Figure 4). The ctgl026 (Figure 5B) is similar 
to a carrot PR which has been related to early embryo 
development, being expressed in the endosperm and 
secreted in the apoplast, thus positively regulating 
embryo fate and patterning [17]. It is interesting to note 
that in cv Springcrest, the down-regulation of ctgl026 at 
S3 and S4 occurs at a slower rate than in Fantasia and 
sir, thus confirming the uncoupling of seed and meso- 
carp development also at the molecular level. The differ- 
ent kinetics observed for the expression of S3 marker, a 
gene encoding a prunin storage protein (ctgl540, Figure 
5C) in sir indicates that in this selection, as well as the 
blocked development of the mesocarp, some variations 
in seed storage accumulation may also exist. The appar- 
ent delay in transcript accumulation measured in 
Springcrest is again due, as in the case of ctg3431, to 
the uncoupled development of seed and pericarp. 
Ctg3563, encoding a LEA (late embryogenesis abundant) 
protein, is a very reliable marker of S4, in both Fantasia 
and sir, indicating that the seed can reach a fully 
matured stage in both genotypes. The very low levels of 
LEA gene expression detected at S4 in Springcrest are 
consistent with the uncoupling that exists between seed 
and pericarp maturation in this genotype. 

Concerning the mesocarp, ctg2909, marking SI and 
S2I, encodes a putative RD22-like protein, whose 
expression in both Arabidopsis and grape is partially 
under the control of ABA and claimed to be involved in 
stress responses [18,19]. Since the levels of this hormone 
in peach mesocarp were shown to follow a biphasic pat- 
tern with two peaks at S2I and S4 [20], the increasing 
expression of ctg2909 at early mesocarp development 
might be related to the level of ABA. However, while 
the hormone also peaks at S4, the expression levels of 
this gene did not, thus indicating a dual regulatory 
mechanism triggering its expression, possibly also under 
a developmental control as shown in the seed of Arabi- 
dopsis [19]. The delayed decay of ctg2909 expression 
observed in Springcrest might be related to the higher 
growth potential of this early ripening variety documen- 
ted by the S2 phase length, which is significantly 
reduced compared to cv Fantasia (Figure 5E). The S2 
phase is marked by ctgl751 (Figure 5F), coding for a 
serine carboxypeptidase (SCP). SCPs are members of the 
a/p hydrolase family of proteins, claimed to function 
also as acyltransferases and lyases in the biosynthesis of 
secondary metabolites [21]. Taking into account that the 
most important event occurring at S2II is endocarp lig- 
nification, an indirect role in this process might be 
hypothesized for ctgl751. Ctgl823 (Figure 5G) was 
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shown to be a good S3 marker in Fantasia, but not in sir 
and Springcrest. Since this gene encodes a putative 
senescence-associated protein, a likely failure of the 
senescence process and/or of its entry phase may be 
hypothesized in the two genotypes in which mesocarp 
development is either slowed down or accelerated. Inter- 
estingly, when mesocarp development is slowed down, 
as in sir, the peak of expression of ctgl823 is antici- 
pated, whereas in the other case (i.e. in Springcrest), in 
which mesocarp ripens very rapidly, the peak is almost 
absent. It may be speculated that an overly precocious 
start of senescence would not allow the fruit to shift 
from maturation to ripening [22], and, vice versa, an 
acceleration of fruit ripening is achieved if senescence is 
not initiated. For the S4 phase, a very reliable marker is 
represented by ctg57 (Figure 5H), coding for an already 
partially characterized peach Aux/IAA protein [10]. Its 
expression was shown to increase at early S4, most 
likely under a developmental control, thereafter decreas- 
ing when ethylene climacteric is fully installed. Accord- 
ingly, ethylene treatments were shown to reduce the 
specific transcripts. Besides fully agreeing with previous 
data, the expression profiles shown here may also repre- 
sent correlative evidence for a putative functional role. 
Indeed, no rise of expression was measured in the meso- 
carp of sir, consistent with the block/slowdown of devel- 
opment and ripening. Moreover, in the case of 
Springcrest, a high ethylene-producing variety [23], the 
rise in expression of ctg57 is both anticipated, parallel- 
ing ripening kinetics, and less pronounced than in Fan- 
tasia, in agreement with a negative effect exerted by 
higher levels of ethylene. 

Possible mechanisms involved in seed-pericarp cross- 
talk should take into account the vascular and cellular 
connections existing between the two organs. It has 
been shown that all the maternal tissues of pericarp and 
seed (integuments) are intensively interconnected (Viz- 
zotto, personal communication), while nucellar tissue is 
excluded from the plasmodesmata network. This implies 
that the flux of metabolites, as well as signaling mole- 
cules between embryo and fruit, must occur through the 
apoplast. Taking into account that hormones play a 
pivotal role in the regulation of seed and fruit develop- 
ment, it has been assumed that they might also be 
involved in the cross-talk between the two organs. The 
heat map data (Figure 6) will therefore be discussed tak- 
ing into account the consistency of the colors in the fol- 
lowing main two-by-two comparisons: ST/H, SRG/H, 
TFs/H, and sTFs/SRG. More specifically, considering 
the first one (ST/H), consistency of colors may indicate 
a relationship between the hormone-related response 
and activation of the corresponding signal transduction 
pathway. In the second comparison (SRG/H), the same 
parameter may provide information about the hormone 



specificity of the transcriptional response, and, at the 
same time, of the possible cross-talk between hormones. 
A double comparison (TFs/H, and sTFs/SRG) may allow 
it to be pointed out if other players besides the TFs are 
involved in the regulation of the downstream processes, 
and if a specific response is mediated by hormone-speci- 
fic TFs. Auxin, cytokinins, and gibberellins are generally 
considered to be the most relevant hormones for early 
seed and fruit development, whereas abscisic acid and 
ethylene play important roles in seed maturation and 
fruit ripening. From the point of view of the cross-talk 
between seed and mesocarp, comparisons should refer 
to the same developmental stage, i.e. ES/EM and LS/ 
LM. Concerning auxin, the data presented here point 
out that the specificity of the response to the hormone 
is higher in ES and LM, although the relationship 
between the overall HORMONOMETER (H) and ST 
data indicates that mesocarp is always more sensitive 
than the seed to the hormone. Taking into account that 
the presence of a viable seed is required for fruit set and 
development in peach [2], and that the overexpression 
of auxin biosynthetic genes in the ovary stimulates the 
parthenocarpic fruit development in several species [6], 
it may be hypothesized that the signal produced by the 
developing seed might be either the auxin itself exported 
to the fruit, as demonstrated in tomato [24], or a sec- 
ondary messenger whose target at the fruit level 
includes a large subset of auxin-responsive genes. This 
is consistent with both the high specificity of the auxin 
response shown here in the early developing seed and 
the higher sensitivity to the hormone displayed by the 
mesocarp paralleled by a strong hormone response. 
Among the mesocarp auxin responsive genes, several 
encode elements regulating transport (ctg2448, ctg2449 
and ctg2789 [25] Additional file 1), indicating that auxin 
movement in this tissue is a relevant process, thus 
strengthening the hypothesis that auxin produced by the 
seed may behave as a signal efficiently transported to 
and within the mesocarp. An Aux/IAA-encoding gene 
(ctg358) showed an opposite transcription profile in the 
two organs, being abundant in ES and LM. It has been 
demonstrated that its tomato orthologue (i.e. LS-IAA9, 
[26]) acts as a repressor of auxin signaling. Thus, its 
expression in young organs (low in mesocarp, high in 
seed) seems to confirm that the hormonal response is 
not at the synthesis site. Finally, the expression of 
ctg2655, a SAUR-like IAA responsive gene [27] was 
found to be higher in mesocarp than in seed (see also 
Figure 3), thus suggesting a higher auxin level in EM 
than in ES [28,29]. 

The main process regulated by CKs is cell division, 
occurring at early development in both seed (endo- 
sperm) and mesocarp. In the former, there is an up-reg- 
ulation of signal transduction elements, such as ctg2370 
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coding for a histidine-containing phosphotransfer pro- 
tein [30] whose transcription is abundant in very young 
organs (Figure 3). The corresponding substantial activa- 
tion of hormonal targets, including several CK-specific 
genes, might differ in the two organs. For example, a 
cellulose synthase (ctg3673) is activated in EM but not 
in ES, cytokinesis being an LS event, whereas cyclin D3 
(ctg779) was up-regulated in both organs at the early 
stage. However, this transcriptional response did not 
just involve CK-specific TFs, implying that other regula- 
tory elements may determine the hormone-specific gene 
activation. A similar activation of signal transduction 
elements to that found in the seed is present in the 
mesocarp at early development. However, the overall 
and the CK-specific target activation are not correlated 
to the hormone action, suggesting that CKs may regu- 
late mesocarp cell division at the post-transcriptional 
level [31], either alone or in cooperation with other phy- 
tohormones. Moreover, considering the inter-organ 
comparison, it is noteworthy that during early develop- 
ment the seed displayed a higher sensitivity to CKs than 
the mesocarp but a lower specificity of response. The 
amount of the overall transcriptional response observed 
in the seed may be due to the involvement of other hor- 
mones besides CKs [32]. During late development, an 
inverse situation was observed compared to the early 
phases. In fact, the high activation of signal transduction 
pathways occurring in the seed was uncoupled from the 
overall transcriptional response, which was even more 
specific in the mesocarp. The CK-mediated up-regula- 
tion of genes encoding sorbitol dehydrogenases (ctg636 
and ctgl378, Additional file 1) appears particularly inter- 
esting, as this might increase the sink strength of the 
seed and attract photoassimilates to the entire fruit, 
which become more competitive in the partitioning pro- 
cess [33]. 

From a physiological point of view, GAs play a stimu- 
latory role in fruit development, as shown by the ability 
to induce parthenocarpy in several species [34] when 
applied in post-bloom phase and/or early development. 
The initial phases of endosperm and embryo develop- 
ment are usually related to a high level of GAs [35], 
while seed maturation is paralleled by a decay of free 
GAs and increase of their conjugates. The HORMON- 
OMETER data confirmed these results both in seed and 
mesocarp, except for TFs in the latter. In fact, the most 
relevant transcriptional response occurred during early 
development at seed level as pointed out by the ES-spe- 
cific expression of ctg3431 (Figure 5) encoding an 
orthologue of the Arabidopsis LTP1 (AT2G34580), 
which is classified as a GA-responsive gene (see at 
http:// genome.weizmann.ac.il/hormonometer/) involved 
in embryo patterning [36]. In the mesocarp, a low corre- 
lation was observed between the TF-related 



transcriptional response and GA action, implying the 
activation of complex regulatory mechanisms that may 
play relevant roles in the cross-talk between seed and 
mesocarp. A possible mode of interaction might be the 
EM specific expression of a gene coding for a Zinc fin- 
ger protein (ctgl87), whose Arabidopsis orthologue 
(AT2G04240, XERICO) interacts with DELLA proteins, 
is repressed by GA, and causes ABA accumulation when 
over-expressed [37]. However, since this transcriptional 
response lacked specificity, it might be hypothesized 
that GA action also depends on the interaction with 
other hormones. It has recently been demonstrated that 
auxin induced parthenocarpy via GAs in unpollinated 
tomato ovaries [38]. Furthermore, the peculiar expres- 
sion profile of ctgl391, encoding a GAST-like protein, 
orthologue of Arabidopsis GASA6 (AT1G74670), in EM 
is confined to S2 and S4 stages, when cell enlargement 
is slow (Figure 3). These data are in agreement with the 
observed inhibition of cell elongation conferred on both 
Arabidopsis seedling and strawberry fruit over-expres- 
sing the Fragaria orthologue FaGAST [39]. During late 
development, in spite of the slight correlation existing in 
terms of GA-specific response, the other gene sets 
appeared not to be correlated to the hormone action. It 
may be deduced from this that the role of GA in the 
cross-talk between seed and mesocarp is negligible dur- 
ing late development. 

ABA is known to play an antagonistic role with 
respect to auxin, GAs and CKs, as observed during fruit 
development in avocado [40] and tomato [41,42]. 
According to the HORMONOMETER, this antagonism 
was largely confirmed in the seed, the transcriptional 
response being correlated with higher levels of the hor- 
mone in LS compared to ES, also when the ABA-speci- 
fic subset was considered. In fact, during late seed 
development, ABA levels are known to increase and 
GA-related genes such as ctg3430, encoding a LTP-like, 
are down-regulated (Figure 3 and Additional file 1). 
This physiological parameter is paralleled by a consis- 
tent transcriptional response in which TFs belonging to 
WRKY (ctgl545), HD (ctg499), Aux/IAA (ctg768), bZIP 
(ctg 724) and DREB-like AP2 (ctg 4674) families are 
involved. Given this interpretation and taking into 
account that during both early and late development 
ABA ST pathways and ABA-target responses are more 
active in the mesocarp, the hormone may play a more 
relevant role in the development of each organ, rather 
than in seed-mesocarp cross-talk. In this context, the 
ABA pool of maternal and zygotic origin may trigger 
independent transduction pathways. 

The well-known role of ethylene in peach ripening 
[9,10] was confirmed by the higher level of transcription 
of its ST elements (ctg4109, ctg244 and 4757 coding for 
an ETR2-like ethylene receptor and two ERFs, 
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respectively) measured at the mesocarp level during late 
development. It is worth noting that ethylene-related 
transcriptional response in the LM/EM comparison 
resembles that of ABA, most likely because of the signif- 
icant number of transcriptional targets shared by the 
two hormones (50 out of 216 ABA- and 235 ethylene- 
responsive genes). This was not observed in the inter- 
organ assessments, in which, differently from ABA, a 
weaker overall transcriptional response was pointed out, 
with the only exception being ethylene-specific TFs that 
were more represented in mesocarp. As regards cross- 
talk, the role of ethylene might be limited to the very 
early phase of fruit development, as demonstrated in 
tomato [43]. This might be a consequence of the fact 
that ethylene acts mainly within the cell where it is 
synthesized. Also in this case, the hormone of maternal 
and zygotic origin may activate independent pathways 
controlling different processes. 

Conclusions 

In this research, genes were identified marking different 
developmental phases of seed and mesocarp. The reliabil- 
ity of these molecular markers was tested in two subse- 
quent years and a further functional validation was 
carried out in three different genotypes. In the latter case, 
data indicate that, while seed markers represent reliable 
tools in all the tested varieties, in the case of the meso- 
carp the different developmental and ripening traits of 
the various genotypes somewhat affect the expression of 
marker genes, consistently with their putative functions 
and cv characteristics. The most critical phases, from the 
point of view of mesocarp marker retrieval, were S2II 
and S3. This might be related to the high divergence in 
pericarp development among the different genotypes, as 
pointed out above. However, this limitation may be par- 
tially overcome by using mesocarp markers as a whole, 
therefore increasing their discriminating power. 

As regards the cross-talk between seed and pericarp, 
possible candidate signals were identified among hor- 
mones. In the early phases, when the cross-talk is more 
vital for fruit set, the candidates are auxin, CKs, and GAs, 
acting either directly (auxin) or indirectly as signals, 
whereas ABA and ethylene appear to be involved later on. 

Further investigations relying upon the availability of 
whole genome platforms will allow enrichment of the 
marker genes repertoire and elucidation of the cross- 
talk mechanisms between the two organs, taking into 
account, besides hormones, other players such as hor- 
mone peptides and microRNAs. 

Methods 

Plant materials 

Fruit growth analysis was conducted on peach trees of 
cv Fantasia grown on the experimental farm of the 



University of Padova (Legnaro), Italy, as described by 
Tonutti et al., 1997. Fruits from 10 trees were collected 
at 42, 60, 81, 106 and 123 days after full bloom (DAFB), 
corresponding to the first exponential growth phase 
(SI), the onset (S2I) and end (S2II) of pit hardening, 
second exponential growth phase (S3) and ripening (S4), 
respectively. For each sampling, mesocarp and seed 
were excised from 30 fruit, pooled in three biological 
replicates and then immediately frozen in liquid nitro- 
gen and stored at -80°C until use. To monitor seed 
development, seeds were excised from fruit at weekly 
intervals from late SI to ripening. Seed and embryo 
length were measured by stereomicroscopy [44]. 

Seed and fruit of two additional genotypes (cv 
Springcrest and selection slow ripening - sir-) charac- 
terized by uncoupling of seed and fruit development, 
were used for the validation of marker gene functions. 
In Springcrest, an early ripening cultivar, fruit ripening 
occurs when seed development is still in progress. In 
fact, seeds become viable only after in vitro cultivation. 
sir is a selection obtained from a free-pollinated popu- 
lation of Fantasia, characterized by a block of meso- 
carp development at stage S3 but with a fully 
developed seed. Sampling of mesocarp and seed was 
performed throughout fruit development as previously 
described. 

Fruit growth analyses were performed in 2008 and 
repeated in 2010; the array experiments were performed 
on samples collected in the former year, whereas expres- 
sion data were validated by qRT-PCR on samples of 
both years. Only data related to 2008 are presented and 
discussed in the paper. Data from 2010 are given in the 
supplementary material (Additional file 3). 

Transcriptome analysis 

For each sampling data, total RNAs were extracted, as 
described in [45], from each of the three biological repli- 
cates of seed and mesocarp and stored at -80°C for tran- 
scriptional analysis. 

To elucidate the interactions between seed and meso- 
carp, a mass gene approach was followed by using the 
uPEACHl.O as described in [10]. Comparisons were 
made by pooling stage 1 and 2 (named early develop- 
ment, E), and stage 3 and 4 (named late development, 
L), separately for mesocarp (M) and seed (S), and using 
a simple loop experimental design (Figure 1). 

Data were analyzed using the TM4 software platform 
[46] as previously described [10]. A SAM (Significance 
Analysis of Microarrays [47]) analysis was performed to 
identify significantly differentially expressed genes using 
a False Discovery rate of 0% (90 th %ile). Among these, up 
and down regulated genes were identified assuming a 
threshold ratio of expression as log 2 higher than 1 and 
lower than -1, respectively. 
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To improve the annotation of targets spotted on the 
uPEACHl.O platform, all the oligo sequences were 
blasted against a transcript dataset obtained by assem- 
bling 280,000 454 reads (Additional file 7) with the 
about 90,000 Sanger Prunus persica ESTs present in the 
NCBI database. The 454 reads have been obtained from 
a normalized library constructed by pooling equal 
amounts of mesocarp RNA from stages SI, S2, S3 and 
S4 (GenXPro GmbH, Germany). The 32,162 new con- 
tigs present in this new database have been compared to 
those used to develop the uPEACHl.O platform and, if 
longer than the old ones, used for BLAST analysis. Con- 
tigs (Additional file 8) were analyzed by BLAST against 
already classified proteins from Arabidopsis (TAIR 10 
release) to categorize them by using the GO terms 
developed by TAIR http://www.arabidopsis.org/portals/ 
genAnnotation/functional annotation/ ontologies.jsp for 
the biological processes ontology. Based on the best 
BLAST search results and using a cut-off e value of l*e- 
10 , the peach genes were assigned to the categories 
according to the most similar Arabidopsis genes (Addi- 
tional file 1). 

Differentially expressed genes were visualized with 
Venn diagrams drawn with Venny [48], clustered 
according to their expression profiles by using the Qual- 
ity Threshold Clustering (QTC) coexpression algorithm 
[49] and grouped in four main charts allowing intra and 
inter-organ as well as developmental-stage comparisons. 

The data discussed in this paper have been deposited 
in NCBI's Gene Expression Omnibus and are accessible 
through GEO Series accession number GSE22582 
http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi? 
acc=GSE22582. 

qRT-PCR was performed and the obtained data 
manipulated as previously described [10]. Briefly, 3 ug of 
total RNA for each sample, pre-treated with 1.5 units of 
DNasel, was converted to cDNAs by means of the 
"High Capacity cDNA Archive Kit" (Applied Biosys- 
tems), which uses random examers as primers. Primer 
sequences for the selected genes are listed in Additional 
file 6. Oligonucleotides PpNlfor (CCAGGAGAATC 
GGTGAGCAGAAAA) and PpNlrev (TCGAGGGTG- 
GAGGACTTGAGAATG) annealing to the peach puta- 
tive transcript ppa009483 m, orthologous to Arabidopsis 
AT4G34270, were used to amplify the reference gene. 
The peach reference gene was selected starting from 
Arabidopsis homologous genes [50], tested for transcript 
normalization in peach (Tadiello and Trainotti, unpub- 
lished results) and chosen to normalize qRT-PCR data 
because of its superior result compared to the previously 
used Internal Transcribed Spacer of the ribosomal RNA 
[10]. Reactions were performed using 10 uL of the 
"Syber green PCR master mix" (Applied Biosystems), 
with 0.05 pmoles of each primer, in the "7500" 



instrument (Applied Biosystems). The obtained CT 
values were analyzed by means of the "Q-gene" software 
([51]), averaging three independently calculated normal- 
ized expression values for each sample. Expression 
values are given as mean of the normalized expression 
values of the triplicates, calculated according to equation 
2 of the "Q-gene" software ([51]). Differences in expres- 
sion values among probes reflect different quantities of 
target amounts. Numerical values obtained with these 
calculations were transformed into graphics or used to 
build heat maps with MS Excel. 

The hormonometer analysis 

The HORMONOMETER http://genome.weizmann.ac.il/ 
hormonometer/ is a bioinformatic tool for assessing any 
transcriptome response according to the perspective of 
similar events occurring upon hormonal activation [13]. 
A vector-based correlation is calculated by comparing 
the variation of the transcriptome in a query experiment 
with an indexed list of pre-calculated transcriptional 
responses established by published hormone treatments 
in Arabidopsis [52]. Input data, for each gene, consist of 
the fold change calculated by directly dividing the nor- 
malized expression values measured for the two samples 
to be compared, and the respective P-value of its signifi- 
cance. When the variations detected in the query resem- 
ble those of the reference pool related to a certain 
hormone, it is assumed that the same hormone may 
have caused the transcriptional response observed in the 
query. In the HORMONOMETER output data, the 
numeral 1 indicates a complete correlation in terms of 
direction and intensity of the hormone index with the 
queried experiment, 0 indicates no correlation, and -1 
indicates the highest possible anti-correlation for each 
transcript in the index [13]. Given that input data for 
each gene derive from a two-by-two comparison (for 
example, sample A versus sample B), correlation and 
anti-correlation indicate higher levels of the active hor- 
mone that are transduced into a measurable transcrip- 
tional response in either sample A or B, respectively, 
whereas no correlation (the numeral 0) indicates that 
the levels of the hormone are the same in both samples. 

Since HORMONOMETER only accepts Arabidopsis 
data as input along with their corresponding locus name 
and Affymetrix probe IDs, the putative Arabidopsis 
orthologues of peach genes, obtained by using the best 
BLAST hit of the updated uPEACHl.O database against 
TAIR 10, were used as input data with peach expression 
values. In addition to the whole set of peach genes, 
three subsets were submitted to HORMONOMETER: i) 
genes with hormone-specific responsiveness (i.e. that are 
not multiple targets of hormones), ii) hormone-respon- 
sive genes encoding TFs, and iii) genes encoding TFs 
with hormone-specific responsiveness. 
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Additional material 



Additional file 1: Average microarray signal values for all the 
uPEACHI.O probes. Column headings are as follows: A: contig_name; B: 
Operon ID (probe identification of the manufacturer); C: oligo sequence; 
D: Annotation source (sequences used to identify the best TAIR 10 BLAST 
HIT. "Original contig" means that the sequence used comes from the 
older database, while "454" indicates that a new sequence, longer than 
the original one, has been used); E: source ID (database from which 
sequences used for the annotation. Nucleotide sequences are available 
as additional file 7 ("454") and 8 ("Original contigs"); F: TAIR 10 best hit 
(against uPEACHI.O transcripts); G: Sequence description (annotation of 
the TAIR 10 best hit); H: uPEACHI.O original annotation; I to P: mean and 
standard deviation of the mean of normalized signal values in the 
different comparisons (ES/EM, LS/LM, LS/ES, LM/EM); Q and R: tags for 
the identification of genes included in the Venn diagrams (Figure 2); S: 
probes positive to SAM in all 4 comparisons. From column T to AA: 
genes involved in hormone metabolism and signaling for auxin, CKs, 
GAs, ABA, ethylene, jasmonate, salicylate and brassinosteroids are marked. 

Additional file 2: Correlation between microarray and qRT-PCR 
expression values. The qRT-PCR expression values of 29 genes, listed in 
Additional file 6, were plotted against the microarray hybridization 
signals and correlation indexes (Pearson coefficient is reported in the 
inserted rectangles) have been calculated separately for each direct 
comparison (LM/EM, blue diamonds; ES/EM, purple squares; LS/ES, green 
triangles; LS/LM, light blue circles) (panel A). The validation for 21 
randomly selected contigs is shown in the panel B. 

Additional file 3: Validation of seed and mesocarp markers in the 
2010 season. Expression profiles of seed (panel A) and mesocarp (panel 
B) selected markers throughout fruit development. The expression values 
are given as a percentage distribution throughout development of the 
total gene expression (100%). Fisher's Least Significant Difference (LSD) 
was calculated for each gene time series using "agricolae" R package (de 
Mendiburu Felipe, A statistical analysis tool for agricultural research. MS 
thesis. Universidad Nacional de Ingenieria, Lima-Peru. 2009). Panel A: 
ctg3431 LSD = 8.6%, ctgl026 LSD = 1 1.2%, ctgl540 LSD = 15.1%, ctg 
3563 LSD = 7.6%. Panel B: ctg2909 LSD = 4.8%, ctg 1751 LSD = 4.6%, 
ctg 1823 LSD = 5.4%, ctg57 LSD = 3.8% 

Additional file 4: Suitability of uPEACHI.O for the HORMONOMETER 
platform. Table 1: Number and percentage of putative hormone-related 
genes spotted on the uPEACHI.O microarray. Table 2: Number of genes 
representing putative common targets of different pairs of hormones as 
assessed in Arabidopsis [13]. The total number of hormone-related genes 
is given in bold on the diagonal. 

Additional file 5: Heat map showing the relationship between the 
expression of signal transduction and hormone target genes Panel 

A. The heat map was produced by considering the genes involved in 
the signal transduction (ST) for jasmonic acid (JA), salicylic acid (SA), and 
brassinosteroids (BR). HORMONOMETER data were grouped into 
hormone-responsive genes (H), genes with hormone-specific 
responsiveness (SRG), hormone-responsive genes encoding TFs (TFs), and 
genes encoding TFs with hormone-specific responsiveness (sTFs). For 
each hormone, the following comparisons have been analyzed: LS/SE, 
LM/EM, ES/EM and LS/LM. Panel B Color codes for ST genes and 
hormone-responsive genes (HORMONOMETER). For ST, red and green 
represent up- and down-regulation, respectively. In the 
HORMONOMETER, orange (value = 1), white (value = 0), and blue (value 
= -1) indicate a complete correlation, no correlation, or anti-correlation, 
respectively, in terms of direction and intensity of the hormone index 
with the queried experiment [13]. 

Additional file 6: List of the oligonucleotides used for the qRT-PCR 
analyses. The "Contig_nos" refers to the peach contig numbers in the 
database used to prepare the oligo probes of the uPEACHI.O microarray. 

Additional file 7: Sequences of 454 reads Nucleotide sequences of 
454 reads. 

Additional file 8: Sequences of Contigs Nucleotide sequences of 
contigs. 
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